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ABSTRACT 

Exotic dark matter together with the vacuum energy (associated with the 
cosmological constant) seem to dominate the Universe. Thus its direct detection 
is central to particle physics and cosmology. Supersymmetry provides a natural 
dark matter candidate, the lightest supersymmetric particle (LSP). One essen- 
tial ingredient in obtaining the direct detection rates is the density and velocity 
distribution of the LSP. The detection rate is proportional to this density in our 
vicinity. Furthermore, since this rate is expected to be very low, one should 
explore the two characteristic signatures of the process, namely the modulation 
effect, i.e. the dependence of the event rate on the Earth's motion and the 
correlation of the directional rate with the motion of the sun. Both of these 
crucially depend on the LSP velocity distribution. In the present paper we study 
simultaneously density profiles and velocity distributions based on the Eddington 
theory. 

Subject headings: Cosmology: Eddington theory, velocity profiles, rotational curves- 
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1. Introduction 

In recent years the consideration of exotic dark matter has become necessary in order 
to close the Universe (Jungman et ah 1996). The COBE data (Smoot et al. 1992) suggest 
that CDM (Cold Dark Matter) component is at least is at least 60% (Gawser et al. 1988; 
Gross et al. 1998) of the total mass. On the other hand evidence from two different teams, 
the High-z Supernova Search Team (Riess et al. 1998) and the Supernova Gosmology Project 
(Somerville et al. 2002) ' (Perlmutter et al 1999; Perlmutter et al. 1997) suggests that the 
Universe may be dominated by the cosmological constant A. Thus the situation can be 
adequately described by a baryonic component Qb = 0.1 along with the exotic components 
^CDM — 0.3 and Ha — 0.6. In a more detailed KCDM analysis (Primack 2002) one finds: 

= 0.040 ± 0.002, = QcDM = 0.33 ± 0.035, ^hdm < 0.05, Qa = 0.73 ± 0.08 

(see also (Turner 1990) and (Einasto 2001)). Since the non exotic component cannot exceed 
40% of the GDM (Jungman et al. 1996), (Alcock et al. 1995), there is room for the exotic 
WIMP's (Weakly Interacting Massive Particles). In fact the DAMA experiment (Bernabei 
et al. 1996) has claimed the observation of one signal in direct detection of a WIMP, which 
with better statistics has subsequently been interpreted as a modulation signal (Bernabei 
et al 1998; Bernabei et al. 1999). 

In the most favored scenario of supersymmetry the LSP can be simply described as 
a Major ana fermion, a linear combination of the neutral components of the gauginos and 
Higgsinos (Jungman et al. 1996; Vergados 1990; Gomez and Vergados 2001; Gomez et al. 
2000a,b; Gomez and Vergados 2001; Anrowit and Nath 1995, 1996; Bottino et al. 1997; 
Bednyakov et al. 1994). 

Since this particle is expected to be very massive, > 30GeF, and extremely non 
relativistic with average kinetic energy T < IQQKeV ^ it can be directly detected (Vergados 
1996; Spira et al. 1995; Kosmas and Vergados 1997) mainly via the recoiling of a nucleus 
(A,Z) in elastic scattering. 

In order to compute the event rate one needs the following ingredients: 

1) An effective Lagrangian at the elementary particle (quark) level obtained in the 
framework of supersymmetry as described , e.g., in (Jungman et al. 1996; Anrowit and 
Nath 1995, 1996; Bottino et al. 1997; Bednyakov et al. 1994). 

2) A procedure in going from the quark to the nucleon level, i.e. a quark model for the 
nucleon. The results depend crucially on the content of the nucleon in quarks other than u 
and d. This is particularly true for the scalar couplings as well as the isoscalar axial couphng 
(Drees and Noijiri 1993b,c; Cheng 1988, 1989). 
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3) Compute the relevant nuclear matrix elements (Ressel et al. 1993; Divari et al. 2000) 
using as reliable as possible many body nuclear wave functions. The situation is a bit simpler 
in the case of the scalar coupling, in which case one only needs the nuclear form factor. 

4) The LSP density and velocity distribution. Among other other things, since the de- 
tection rates are expected to be very small, the velocity distribution is crucial in exploiting 
the characteristic experimental signatures provided by the reaction, namely: a) the mod- 
ulation of the event rates due to the earth's revolution around the sun (Vergados 1998, 
1999)^(Vergados 2000) and b) the correlation of the rates with the Sun's direction of motion 
in directional experiments, i.e. experiments in which the direction of the recoiling nucleus 
is observed (Vergados 1990; Buckland et al. 2000). To obtain the right density and velocity 
distributions is the purpose of the present paper. 

In the past various velocity distributions have been considered. The most popular one 
is the isothermal Maxwell- Boltzmann velocity distribution with < >= {3/2)vq where 
vq is the velocity of the sun around the galaxy, i.e. 220 km/s. Extensions of this M- 
B distribution were also considered, in particular those that were axially symmetric with 
enhanced dispersion in the galactocentric direction (Drucker et al. 1986; Vergados 2000). In 
such distributions an upper cutoff Vgsc — 2.84i;o was introduced by hand. 

Non isothermal models have also been considered. Among those one should mention 
the late infall of dark matter into the galaxy, i.e caustic rings (Sikivie 1999(, 1998; Vergados 
2001a; Green 2001; Gelmini and Gondolo 2001), and dark matter orbiting the Sun (Copi et 
al. 1999). 

The correct approach in our view is to consider the Eddington approach (Eddington 
1916), i.e. to obtain both the density and the velocity distribution from a mass distribution, 
which depends both on the velocity and the gravitational potential. This approach has been 
extensively studied by Merritt (Merritt 1985) and recently applied to dark matter by UUio 
and Kamionkowski(Ullio and Kamiokowski 2001) 

2. Density Profiles 

As we have seen in the introduction the matter distribution can be given as follows 

dM = 27r /($(r), Vr, Vt) dx dy dz Vt dvt dvr (1) 

where the function f the distribution function, which depends on r through the potential 
$(r) and the tangential and radial velocities and Vr (we assume axial symmetry in velocity 
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space) . Thus the density of matter p satisfies the equation 

dp = 2tt f{^{r),Vr, Vt) Vt dvt dvr (2) 

It is more convenient instead of the velocities to use the total energy E and the angular 
momentum J via the equations 

J^vtr , 2E^vl + — + 2^r) (3) 

The use of these variables, which are constants of motion, is very useful, when one wants to 
study equilibrium states. We thus find 

2. I I f(E,J)J ^^^^ 

The hmits of integration for E are from $ to and for J from to [2r^(£' — $(r))]^/^. 
Following Eddington we will choose a distribution function of the form 

f{E,J)^K,{-2E)^ (5) 

{E is negative for a bound system), where A is a parameter, which will depend on the type 
of matter, and is a normalization constant, which will be related to the density at some 
point. With this choice of the distribution function it is quite straightforward to find the 
relationship between the density p and the potential. The result is 

p = irA2^+3/27r|$(r)|^+3/2 ^(A + 1,3/2) (6) 

with 

r(a) r(6) 

L [a + 0) 

Eq. (6) can be cast in the simple form: 

P = Pa(0)(x(^))^+^/^ (8) 

with 

, , $(x) , , , 

X{x) = X = r/vs (9) 

with $0 — ^(0) and the galactic radius (position of the sun). The constant is related 
to the density at the origin via the relation: 

7r(2|$o|)^+^/W + 1,3/2) ^'^^ 
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The above distribution function is isotropic. Following Merritt (Merritt 1985) we can intro- 
duce an anisotropy by modifying the distribution function as follows: 

f{E,J)^K,{-2E)\l±:^) (11) 

'a 

Instead of the parameter we find it convenient for our applications later on (see below) 
to adopt the more recent conventions and write the above equation as follows: 

f{E, J) = K^{-2E)'[1 + a.^^^] (12) 

where Vm is the maximum velocity allowed by the potential, to be specified below, and ag 
the asymmetry parameter. Proceeding as above we find that this induces a correction to the 
density of the form: 

Ap = i^A-2^+'/'vr|$(r) 1^+^/2 + l,5/2)^x^ (13) 
3 K 

Combining Eqs (6) and (13) we get 

px{r) = px{0)Mx) (14) 

M^) = ix{x)mi + I a f|^$^^^x(^)] (15) 



with a = as\^o\/vl^. 

Since the scale of the potential appears only via the parameter a one, in principle, could 
have two a parameters, one for Matter (a^) and one for dark matter (adm). In the present 
work wc will assume that they are equal. We remind the reader that ag is the asymmetry 
parameter to be treated phenomenologically. 



3. Allowed Density Functions 

The central question is to specify density ipx{x), the potential x(x) and the mass density 
distribution entering Eq. (1). One can adopt one of two procedures: 

1) Start out with a given density, obtained, e.g., phenomenologically, and find the 
potential by solving Poisson's equation. This way one obtains, at least parametrically, e.g. 
with the radial coordinate as a parameter, the proper relation between the density and 
the potential. This approach has been adopted by a lot of researchers, see e.g. (Widrow 
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2000), (Evans 1994), (Henrriksen and Widrow 1994), (UUio and Kamiokowski 2001). Prom 
this one can obtain the mass density function f{^,v). This approach leads, in general, 
to non analytic, i.e. complicated, velocity distribution, hard to implement in dark matter 
calculations (UUio and Kamiokowski 2001), especially for not spherically symmetric mass 
distributions. 

2) Use the above simple relation between the density and the potential and solve the 
differential equation resulting from Poisson's equation and thus obtain each one of them. 
One hopes that this way one will get a density distribution, which describes adequately 
dark matter distribution globally, via fitting the rotational curves, and in our vicinity. This 
approach can easily deal with asymmetric density distributions. The simplicity of the relation 
between the density and the potential thus yields the bonus that the obtained velocity 
distribution is analytic and can easily implemented in obtaining the event rates for direct 
dark matter detections.lt can also be extended to include more than one power in A. In 
this paper we will follow this approach, and leave it for the future to use semi-analytic mass 
distributions obtained in the approach outlined in the previous paragraph. 

Por a spherically symmetric potential Poisson's equation leads to a differential equation 
of the type: 

xx"{x) + 2x{x) = -A X {x{x)ni + -^x^ix)] (16) 
with n = A + 1/2 and 

^ ^ _A7rGNrjpx{0) ^^^^ 

$0 

with Gjv Newton's gravitational constant. The dimensionless quantity A is assumed to be 
positive (attractive potential). Since the asymmetry parameter is assumed to be small, 
the first term in the right hand side dominates. Introducing the variable ^ = y/Xx and 
x{x) — u{x\/X) we arrive at: 

e«"(0 + 2«'(0 = -e(«(0r[l + ;^eV0] , b = ^ (18) 

The last equation is a non linear differential equation, which must be solved for ^ > 0, with 
the conditions u{0) — 1 and u'{0) — 0. We also demand that the solution remains positive, 
i.e. the solution drops from unity to zero. Por as — the density is a rapidly decreasing 
function of ^. 

The above equation can be solved analytically in the special case a — and only when 
A = — 1/2 (n=l) or X — 7/2 (n=5). In these cases the solutions are: 



(19) 
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and 

A = 7/2^ "(g- d+^VB)'/^ (20) 

The case A = —1/2 may not have physical meaning, but it will serve as an illustrative 
example for the realistic cases to be discussed below. 

We should stress that the last solution associated with A = 7/2, admissible over all 
space, can be considered as providing an adequate description of ordinary matter^. In the 
absence of asymmetry no constraint is imposed on the parameter A by the positivity and 
finiteness condition of the solution. The solution given by Eq. (19) is constrained in the 
range < ^ < tt. If we demand that the range of the potential for dark matter extends far 
out in the halo, e.g. up to Xr — r/vg — 20, we must have A = [tt/x^]^ = 0.025. This gives a 
useful constraint between $o ^-nd px{0). 

The differential equation (18) can, in general, be solved only numerically. For A = 1/2 
we find that there is a singularity at ^ = 12.43, while the first zero of the potential occurs at 
^ = 4.35. This leads to the constraint A = (4.35/20)^ = 0.047. For A = 1 the corresponding 
value is A = 0.065. 

The above results are significantly modified when the asymmetry parameter is turned 
on. One now can see that, depending on the parameters A and as, the shape of the density 
is significantly modified. 

We will examine the special case of A = 7/2, which corresponds to the Eddington 
solution for ordinary matter. As we have seen in the absence of asymmetry the solution can 
be found exactly. In the presence of asymmetry we distinguish two cases: 

1) Positive values of a^. 

For small values of b the solution has no roots. Beyond a critical value of b the solution 
attains the value zero. Then further out for still larger ^ it becomes negative. At some point 
it becomes singular. The extracted values of A range between 1.5 and 0.12. 

2) Negative values of ag. 

In this case when the absolute value of b becomes sufficiently large the solution goes 
through zero. Again a value of A can be extracted. The situation is quite unstable. Generally 



"^The designation of A = 7/2 as rcprc^scinting ordinary matter is due to the work of Eddington (Eddington 
1916), where he points out that this vahio loads to the Plummer distribution. Furthermore this distribution 
provides an adequate description of the galaxy (of ordinary matter) he was concerned with. A discussion of 
some of the thermodynamical aspects of the A = 7/2 distribution is also given in this paper. 
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speaking for arbitrary values of b the solution blows up to infinity without going through 
zero. One cannot extract values of A in this case. 

In any case, these observations are intended as illustrations of what is expected for dark 
matter where the exact solution cannot be obtained. In the case of ordinary matter we will 
not consider ctg different from zero, since the range of the solution is small ( the density is 
assumed to vanish for r > r^). 

Let us return to Eq. (12) and express it in terms of the velocity in our vicinity. We 
obtain: 

f{v) = 7V[-2$(r,) - v']\l + (21) 



where is a normalization constant, which depends on A, and Vm- From this wc see that 
the maximum velocity in our vicinity , Vm, depends on the value of the potential, i.e. 

vl, = -2 $(r,) = -2 % uiVX) (22) 

This means that 

b^^ = (23) 

A 2Am(VA) 

Our strategy is clear: 

1) Vary b so that the solution has its first zero at some value a-nd is monotonically 
decreasing up to that point. 

2) Determine A so that corresponds to a range Xr of the potential xi^)- 

3) From Eq. (23) obtain an acceptable (self-consistent) value of a^, = 2A6m(\/A). 

4) The parameter Vm is obtained via the relation: 

,«(v^) p,(0) 



vl= = -2 1>„u(V-X)^8.G,r r^ (24) 



or equivalently 



'^m — '^o\ffhhy{n,a,K) , hy{n,a,h) — \ ^ i^-V^'^ , n = A-|-3/2 

"^'"^ ^ ^ -V , , ; ^A((^x/A))«[l + |^«(x/A)]^ ' 

(25) 

where vq is our velocity of rotation around the center of our galaxy (220 km/s). In obtaining 
the last equation we used our solution to relate the density Pa(0) to the density of dark matter 
in our vicinity, poi via the equation: 

Pa(0) = 7= ^ , n = A + 3/2 (26) 
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The parameter rj^ is defined by rj^^ — po/(0.3 GeV/cm^). r]^ is normally taken to be unity 
(Jungman et al. 1996). 

It is clear that in our treatment for a given A we have three independent parameters, b, 
Xr and po- The obtained results for some interesting typical cases of ctg are shown in Table 
1. The value of the parameter has the meaning of yesc (see section 6) and it is rewarding 
that it is not far from the value of 2.84 deduced phenomenologically (Druckcr ct al. 1986). 
The potentials obtained for the various such parameters are shown in Figs 1-4. It is clear 
that the shapes of these potentials are similar. 

The densities obtained with the same set of parameters are shown in Figs 5-8. We 
observe that, even though the potentials are similar, the densities differ substantially. The 
asymmetry parameter has a big effect on the density, especially for greater than zero, 
when it changes substantially around x = S. For negative there is small reduction of the 
density in the region of interest, but, as we have already mentioned, only special values of h 
are acceptable. A detailed rigorous mathematical study of the behavior of the solutions for 
h negative is currently under study and it will appear elsewhere. 



4. Rotational Velocities 

We are now in position to obtain the rotational velocity curves. The rotational velocity 
is given by 

vx{x) = (27) 

with ^ 

9x{x) r {x'f dx M^') (28) 

X Jo 

We see that the scale of the rotational velocity is set by ^ in agreement with the data, see, 
e.g., the recent review (Jungman et al. 1996). For dark matter the above integrals can only 
be performed numerically. The obtained rotational velocities are shown in Figs 9-11. For 
ordinary matter (A = 7/2) we make the choice of A = 3 to get the famihar solution 

Xi^)- (1 + 12)1/2 ' 4^7/2 ^IPix)^ (1+^2)5/2 (29) 

(We will drop the index 7/2 and use the index m, for matter, when necessary). Thus the 
rotational velocity due to ordinary matter is now given by 

Vmix) = ^^/gm{x) (30) 
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2.5 5 7.5 10 12.5 15 17.5 20 



x-> 



Fig. 1. — The potential function in units of $o ior A = 1/2. The notation of the curves is: 
Thick sohd hne,intermediate thickness sohd hne, dotted hne, fine dashing, long-fine dashing, 
intermediate-fine dashing and long-long dashing in the order of the asymmetry of Table 1 
(here = 0.000, -0.021, -0.038, -0.080, 0.062, 0.100, 0.150, 0.200 and similarly for the other 
values of A). 
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Fig. 2. — The potential function in units of $o for A = 1. 




Fig. 3. — The potential function in units of $o for A = 



-1/2. 
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Fig. 4. — The potential function in units of $o for ordinary matter. 
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Fig. 5. — The density in units of p(0) for A = 1/2. 
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Fig. 6. — The density in units of p(0) for A = 1. 




Fig. 7. — The density in units of p(0) for A = 



-1/2. 




Fig. 8. — The density in units of p(0) for ordinary matter, \ — 7/2. 
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Fig. 9. — The rotational velocities in units of ^ for A = 1/2. 
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Fig. 10. — The rotational velocities in units of ^ for A = 1. 
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Fig. 11. — The rotational velocities in units of ^ for A = —1/2. 
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For X < 1 we write 
while for x > 1 we have: 

We now distinguish two cases: 
i) Spherical galaxy. 
Using the above density we find: 



9m{x) = g{x) (31) 



= - ^(1) (32) 

X 



^^^^ = l {l + x^f/^ ^^^^ 



ii) Spiral galaxy (disk). 



In this case we will use the same density as above. In the Eddington theory the equation, 
which relates the potential and the density is no longer of the above simple form. One finds: 

to 12' V(0)' ' l + ip/po)" " 

This complexity is of no concern to us since wc do not need to use this complicated equation. 
In the case of ordinary matter we do not need to obtain the velocity distribution. Thus with 
the assumed density we find: 

The above functions are plotted in Fig. 12. In the graphs we show not only the physically 
interesting case, in which the density vanishes outside the radius of the galaxy, but, for 
illustrative purposes, the case of infinite extent of the same density function . We see that, 
since the density falls very fast, it makes very little difference, which form we use. 

By comparing the above graphs of dark matter with those of ordinary matter we see that 
in the Eddington theory ordinary matter can make a significant contribution to the rotational 
velocities, only if the density at the origin is much bigger compared to that of dark matter. 
As we have mentioned above for ordinary matter there is no constraint between $o and po- 
So we simply rescale ordinary matter by a factor Cmdm and write: 



9matter{x) — Cmdm9m{x) (36) 

the scaling factor can be determined by a comparison to the experimental rotational curves. 
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Fig. 12. — The rotational velocities in units of ^ ordinary matter. Thick lines are associated 
with a spherical galaxy, while the thin lines are associated with a spiral galaxy. The finite 
extent of the matter density gives essentially the same results with that of infinite range. 
This happens, because the density falls relatively fast. 
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Thus the rotational velocity due to both matter and dark matter is given by 

Vx{x) = ^^gx{x) + gmatter{x) (37) 

The obtained rotational velocities for two values of Cmdm are shown in Figs 13 - 16 The 
larger value of Cmdm seems to be in better agreement with the data, see, e.g., the recent 
review (Jungman et al. 1996). 

It is clear from these curves that the presence of asymmetry has a dramatic effect on 
the rotational velocities. Thus, in the context of the Eddington theory as discussed here, a 
large asymmetry is excluded by the data on rotational velocities. A detailed fitting of our 
input parameters to the rotational curves will not be attempted here. 



5. Velocity distribution with respect to the galactic center 

The above density via the Eddington formula leads to a velocity distribution of the 
form: 

F(t;,r) oc [-2$(r) -i;Y(l + a,^) (38) 

The above velocities and the distance r are defined with respect to the center of the matter 
distribution, i.e to the center of our galaxy. We note in the context of the Eddington theory 
the velocity distribution cannot be Maxwellian. For a given distance it goes to zero at the 
boundaries of the corresponding ellipsoid. It is customary to consider the value of the above 
distribution in our vicinity, r — Vg- This way it reduces to the product of the local density 
and the velocity distribution. The latter is given by Eq. (21) where is a normahzation 
constant, which depends on A, ccg and Vm- The above notation was introduced to make 
the last equation coincide with the standard expression when the function f is chosen to be 
MaxweUian, i.e. 

^-p(-^)-(l + «4) Exp{-^^) 

for sufficiently small a^. In this limit we see that coincides with the parameter —A 
of Vergados (Vergados 2000) et al Drucker et al. (1986) (in the present work A is used for 
another purpose). 

It is straightforward to find that the normalization factor N is given by 

N-\\ as, Vm) = 27rvl'^'f3{X + 1, 3/2)[l + ^ a, ||^^{^)] (39) 
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Fig. 13. — The rotational velocities in units of for A = 1/2 in the case of spiral galaxy 

with Cmdm = 30. 
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Fig. 14.— 



The same as in Fig. 13 for Cmdm — 2.5. 




Fig. 15. — The same as in Fig. 13 for a spherical galaxy. 
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Fig. 16. — The same as in Fig. 14 for a spherical galaxy. 
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with Vm given by Eq. 22. In the special case of dark matter (A = 1/2) it becomes 

N-\l/2, a„ vj = + ^a,] (40) 

From the above formulas wc see that the velocity of dark matter with respect to the galactic 
center ranges from to a maximum speed Vm = |2$(rs)|^/^. Since the escape velocity is 
determined by the potential $(rs) due to all kinds of matter, the velocity Vm is not simply 
related to Vgsc- Thus, since the distribution function must remain positive, if < its 
absolute value is bounded. This imposes a constraint, since in the traditional analysis with 
only axially symmetric Gaussian distribution it leads to negative a^, i.e. enhanced dispersion 
in the galactocentric direction, a phenomenologically preferred result (Drucker et al. 1986). 
The data on rotational curves may provide an additional constraint on the negative values 
of as- For positive values of ag the constraint coming from the rotational curves is, as we 
have seen, more stringent. 

From then on one proceeds in the usual way to obtain the velocity distribution with 
respect to the laboratory. 

6. Velocity distribution with respect to the laboratory 

For this transformation one needs the velocity of the sun around the galaxy vq — 
220Km/s, a fraction of the escape velocity, which is Vesc — Q25Km/s — 2.84 vq (Drucker et 
al. 1986). 

It is convenient to choose coordinate system with its polar z-axis in the the direction 
of the disc's rotation, i.e. in the direction of the motion of the the sun, the x-axis is in the 
outward radial direction and the y-axis is perpendicular to the plane of the galaxy. Since 
the axis of the ecliptic (Kosmas and Vergados 1997). lies very close to the y,z plane the 
velocity of the earth around the the center of the galaxy is given by 

Ve — f oz + Vi — vqZ + ( sin q; X — cos a cos 7 y + cos a sin 7 z ) (41) 

where a is the phase of the earth's orbital motion, a = 27r(t — ti)/TE, where t is the time of 
observation and ti is around second of June and = 1 year. The magnitude of the Earth's 
velocity is much smaller than that of the sun, i.e. Si = 2vi/vq = 0.27 The velocity of the 
earth around its own axis is even smaller and it is usually neglected. 

One can now express the above distribution in the laboratory frame by writing v' = 
V + Ve , where the prime indicates the velocity with respect to the center of the galaxy. 
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Indicating by y the velocity of the LSP (Lightest Supersymmetric Particle) in units of vq, 
i.e. by defining y — v/vq, we find 

f{y,9,cj>) = N{X,as,ym) [vl - Y{y,eA)]\l + as{Y{y,e)-{ysine cos0-| sina)^)] (42) 
with 

N{\ a„ ym) - vl^+^N{\ a„ v^) (43) 
where iV(A, a^, Vm) is given by Eq. 39, ym = Vm/vQ and 

^2 

Y{y,9,(j)) = 1 + ^ + y'^ + 2y cos9 + 6i[y cos9 cos a sin 7 

— y sin 6* sin cos a cos j + y sin 9 cos sin a] (44) 

with In the conventional axially symmetric Gaussian velocity distribution this function is 
given by 

/(t/,^,0) = ^(«^^^e3c) Exp[-{as + l)Y{y,9,cf>) (45) 
+ as{y sinO cos0 — — sina)^)] 

In this case < y < yesc, but the upper cutoff is introduced here artificially. The normal- 
ization here is defined so that N(as — 0, yesc — 00) = 1 (Vergados 2000) 

The detection rate in direct dark matter experiments is obtained by convoluting the the 
relevant cross section with the above velocity distribution. If the dark matter candidate is 
the LSP, the a-dependence of the above distribution , present only when Si ^ 0, gives rise 
to the modulation effect, i.e. the dependence of the rate on the Earth's motion. This signal 
can be used to discriminate against background. 



6.1. The non directional rate 

The non directional differential event rate is given by: 



dR 

du 



R^T{u) , Tin) = a'\F(u)\'^'(aV^) (46) 



for the coherent mode and 

dR 
du 



= R^Tspin{u) , Tspin = a'^\Fuu)\-^{a^/u) (47) 
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where R and Rgpin are the total rates for the coherent and the spin contributions associated 
with an average LSP velocity, y/ < v'^ > — -sJJ^jTivQ. These parameters, which carry the 
dependence on the SUSY parameters, are the most important ones, but they are not of 
interest in our present calculation. F{u) is the form factor, entering the coherent scattering 
and Fii{u) is the spin response function entering via the axial current (Vergados 2001b). 
The function ^' depends on the LSP distribution velocity employed and is a function of the 
energy Q transferred to the nucleus 

, go = 4.1 X 10^ A-^/^ KeV (48) 
where A is the nuclear mass number and the parameter a is given by: 



a = [V^iXrhvo]'^ (49) 

where //^ is the reduced mass of the LSP-nucleus system and h is the (harmonic oscillator) 
size parameter. 

The function, which is basic to us, is given by 

*(x)= / dy sm9d9 dcf> y f {y , 9 , <i>) (50) 
Jo Jo Jo 

with < X < ym — i + {Si/2) cos a sin /3. 
The total rate is given by: 

/ "^du (51) 
Ju du 

where Umin is determined from the cutoff energy of the detector and Umax — {Um/o)'^ 

From now on we will specialize our results to the case A = 1/2, which yields an adequate 
description of the rotation curves in the case of dark matter. In the case of this velocity 
distribution, unlike the Gaussian one, one cannot approximate the distribution by a power 
series in 5i. The reason is that there may be threshold problems, when the argument of the 
square root approaches zero. To simplify matters we will still make use of the fact that the 
velocity of the Earth around the Sun is much smaller than the velocity of the Sun around 
the galaxy, 5i « 1. So, if we expand the previous expression into a Fourier series with 
respect to the phase of the Earth, a, only the lowest terms will become important. In other 
words to leading order in 5i it can be put in the form: 



R — R[Rq-\-{Ri cosa sin7 — R2 cosct COS7 + R^ sinQ;)5i/2] (52) 
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It turns out that the expansion coefficients R2 and Rs are zero. We can thus conveniently 
fit the rate with the formula: 

Rt [1 + h cosa] (53) 

where h is the modulation amplitude (the difference between the maximum and the minimum 
is equal to 2| 



In the case of no modulation, 5i — 0, the angular integrals can be done analytically to 
yield: 

^{x) = 2nN{l/2,as,ym /^^'^'^'^^ Mx,ym) - ^Mx,ym)] (54) 
with < X < y^n — 1 and 

Jni^j Vm) — Jinti'lT'j Urm X — 1) — Jintin, yrm X + 1) + 2 Jinti^, yrm 1) (55) 

ry 

Jint{n, ym, y)= [Vm " z^]''^/^dz 
Jo 



(56) 



Jint{n, ym, y) = yiVrnT 2Fi{-, —, -, jy) (57) 



The above integral can be done analytically to yield: 

1 -n 3 y2 
2' ^'2'^ 

where 2-^1 is the usual hypergeometric function. For the cases of interest to us here the 
hypergeometric function can be simplified to yield: 

Jint{3,ym:y) = ^Myi - y^' + 3yiy{yiy' - y')'/' + 3yism-\y / y^)] (58) 

Jint{5,ym,y) = Yi^Hyl-yr^' + 5yly{yl-yY' m 

+ I5yiy{yl - y'f''' + 15y^sin-i(|//|/^)] 

J^nt{^,yra,y) = ^[2Ay{yl - y^ + 2^yly{yl - y^f'^ + 3byl,y{ylf'' (60) 
+ lQbyly{ylf'^ + 105y^sin-^(y/|/„)] 

The corresponding expressions for the Gaussian expressions for as cannot be done analyti- 
cally. For the symmetric case, ctg = 0, one finds that: 

m[x) = herf{x - 1) - erf{x + 1)] + 2 er/(l) (61) 

with erf{x) the error function: 

er/(y) = - re-''dt (62) 

Jo 
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Gpsi (x) 




0.5 1 1.5 2 2.5 

Fig. 17. — The function ^{x) for dark matter in the case of the symmetric Gaussian distri- 
bution. 



Rpsi (x) 




x-> 



Fig. 18. — The function "^(x) in the case of dark matter for the choice A = 1/2 of the 
Eddington theory. The graphs have been labeled as in Fig. 9. 
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The above functions are plotted in Figs 17 and 18. 

Prom Fig. 18 we see that both the value of the function ^{x) as well as its range depend 
crucially on the parameters of the model. In the absence of asymmetry the Gaussian model, 
with an appropriate escape velocity put in by hand, yields results, which are almost four 
times larger than those of the Eddington theory. In the presence of asymmetry our results in 
the Eddington theory are substantially larger. This is due to the larger peak value attained as 
well as the larger allowed range of x. On the other hand wc know that in the Gaussian model 
the total rate is not significantly affected by the asymmetry (Vergados 1999). The reason 
for the strong dependence of the results in the present model on the asymmetry is not, of 
course, the asymmetry per se, but the fact that all parameters change with it. In particular 
the value of ym (see table 1). In the Gaussian model the introduction of asymmetry did not 
affect the velocity distribution in any other way, e.g. it did not affect the cutoff value of 
the velocity distribution. The large upper values of x in the function ^'(x), allowed in the 
present model for large Um, are, of course, somewhat controlled by the nuclear form factor. 
In order to get a better feeling for such an effect on the rate, we also plot the function T[u), 
which is proportional to the differential non directional rate. For illustration purposes we 
have chosen to present results for the popular target ^^"^I and a typical LSP mass of lOOGeV 
(see Fig. 19). The results presented are for the coherent mechanism, but we expect very 
small changes when the spin contribution is considered. It is clear that the introduction of 
asymmetry has a profound effect on the rate. We have seen, however, that the large positive 
values of the ag can, in principle, be eliminated from the data on the rotational curves. 

Effects like those discussed above may be more pronounced in the case of modulation, 
not studied in this work. 



6.2. The directional rate 




(63) 



with 




dy sinedd / d(t>y fiy, 9, 4>) X H{X) 



(64) 



Jo Jo Jo 



with H{X) the well known Heaviside function and X is given by 



X — cos © cos ^ + sin © sm^[sin $ sin + cos $ cos 0] 



(65) 



-34- 



where © and $ describe the direction of observation e 

e = sin 0(cos $63; + sin$ej^) + cosBe^ 

(There should be no confusion of the angle $ used here with the potential $(r) used earlier). 
Note the presence of the factor of l/(27r), since the azimuthal integration of the recoiling 
nucleus is not present and we intend to use the same nucleon cross-section both in the 
directional and the non directional case. 

The total rate is proportional to: 

Rdir = / {-r-)dirdu (66) 



mm 



du ' 



Taking again the lowest Fourier components of the obtained rate as a function of the phase 
of the Earth we get an expression similar to Eq. 52. Thus, to leading order in 5i, we can fit 
the total rate by an expression of the form: 

Rdir = 77- tdir [1 + {hi - h2) cos a + hs sin a] (67) 

ZTT 

The parameters tdir and hi i — 1,2,3 are obviously functions of the direction of observation, 
i.e. © and If one observes in the direction of the Sun's velocity h2 — hs — 0. Similarly if 
one observes in a plane perpendicular to the Sun's velocity hi — 0. Instead of tdir it is best 
to use the ratio: 

^ ^ 27r ^ = ^ (68) 
R t ^ ' 

The parameter k is essentially independent of the LSP mass, the nuclear parameters and 

the asymmetry parameter as- But it depends strongly on the direction of observation and is 

expected to correlate strongly with the angle between e and the Sun's direction of motion. 

This correlation provides a an experimental signature perhaps better the modulation with 

the Earth's motion in non directional experiments. 

For 5i 7^ the above integrals over y, 9, 0, especially in the directional case, can only be 
done numerically. Such results will appear elsewhere (Braun et al. 2002). 



7. Conclusions 

In the present paper we studied the density and velocity distributions of cold dark matter 
in the context of the Eddington theory, considering not only symmetric but axially symmetric 
distributions as well. In our approach we used standard simple distribution functions of the 
energy and angular momentum. This lead us to simple relations between the density p 
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and the potential p and $ were obtained by solving numerically poisson's equation in a 
suitable region of space. This procedure allowed us to determine the maximum permitted 
dark matter velocity. Then we were able to study both the rotational curves as well as the 
velocity distribution of dark matter in our vicinity. We should mention that this distribution 
is not Maxwellian and has an upper velocity cutoff built into into it and not put in by hand 
as in the traditional treatment with Gaussian distribution. 

Our results depend on a minimum set of parameters, which were treated as free (see 
table 1). 

We saw that in the context of this theory the predicted rotational velocities , see Figs 
13 - 16, depending on the input parameters, can vary significantly in the presence of the 
asymmetry. This comes mainly from the factor x"^ in the second term of Eq. 15, when the 
parameter a is reasonably large. By comparing our results to the observed rotational curves 
(Jungman et al. 1996), one may constrain the parameters of the model. The best choice 
seems to be the case with a small asymmetry parameter a^. 

We have also made a preliminary study of the effect of the new velocity distribution on 
the direct detection rates for cold dark matter. For illustration purposes wc have selected 
the case A = 1/2. We have seen that, in the context of the Eddington approach, the total 
rates, unhke the case of the Gaussian distribution (Vergados 1999, 2000), sensitively depend 
on the asymmetry parameter a^. This is due to the fact that, when the asymmetry changes, 
the upper value of the velocity distribution also changes. 

It is thus not surprising that, in the Eddington theory, the total (non directional and 
non modulated) event rates for direct LSP detection maybe substantially different from 
those of the phenomenological Gaussian distributions (compare Figs 17 and 18). The strong 
dependence of the rate on the asymmetry parameter remains even after the nuclear form 
factor has been incorporated (see Fig. 19). Results of more detailed calculations of the event 
rates will appear elsewhere. 

The dependence of the directional and/or modulated rates on the velocity distribution 
is currently under study. We expect this dependence to be more pronounced than on the 
total rates. 

8. acknowledgments 

J.D.V. would like to thank the Physics Department of UNISA and Professor S. Sofianos 
for their hospitality. D.O appreciates the hospitality provided by the University of loannina. 



-36- 



REFERENCES 

Alcock, C. et al. 1995, Phys. Rev. Lett. 74 , 2967. 

Arnowitt, R. and Nath, R 1995, Phys. Rev. Lett. 74, 4952. 

Arnowitt, R. and Nath, P. 1996, Phys. Rev. D 54, 2394; hcp-ph/9902237. 

Bednyakov, V.A., Klapdor-Kleingrothaus, H. V. and Kovalenko, S.G. 1994 Phys. Lett. B 
329, 5. 

Bernabei, R. et al. 1996, Phys. Lett. B 389, 757. 
Bernabei, R. et al. 1998, Phys. Lett. B 424, 195. 
Bernabei, R. et al. 1999, Phys. Lett B 450, 448. 
Bottino, A. et al. 1997, Phys. Lett B 402, 113. 

Braun, M., Owen, D., Sofianos, S. and Vergados, J.D. press, to be pubhshed. 

Buckland, K.N., Lehner, M.J. and Masek, G.E. 2000, in Proc. 3nd Int. Conf. on Dark Matter 
in Astro- and part. Phys., Dark2000, Ed. . Klapdor-Kleingrothaus, H.V., Springer 
Verlag. 

Cheng, T.P. 1988, Phys. Rev. D 38, 2869. 

Cheng, H-Y. 1989, Phys. Lett. B 219, 347. 

Copi, C, Heo, J. and Krauss, L. 1999, Phys. Lett B 461, 43. 

Dawson, S. 1991, Nucl. Phys. B359, 283. 

Divari, P.C., Kosmas, T.S., Vergados, J.D. and Skouras, L.D. 2000, Phys. Rev. C 61, 044612- 
1. 

Djouadi, A. and Drees, M. 2000, Phys. Lett. B bf 484, 183. 
Drees, M. and Nojiri, M.M. 1993, Phys. Rev. D 47, 376. 
Drees, M. and Nojiri, M.M. 1993, Phys. Rev. D 48, 3843. 
Drees, M. and Nojiri, M.M. 1993, Phys. Rev. D 47, 4226. 
Drucker, A., Freeze, A., Spergel, D. 1986, Phys. Rev. D 33, 3495. 



Eddington, A. S. 1916, NRAS 76, 572. 

Jaan Einasto, in Dark Matter in Astro and Particle Physics, p. 3, Ed. H.V. Klapdor- Klein- 
grothaus. Springer- Verlag, Berlin-Heidelberg, 2001. 

Evans, N.W., MNRAS, 267, 333. 

Gawiser, E. and Silk, J. 1988, Science 280, 1405. 

Gelmini, G. and Gondolo P. 2001, Phys. Rev. D 64, 623504. 

Gomez, M.E., Vergados, J.D. 2001, Phys. Lett. B 512, 252; hep-ph/0012020. 

Gomez, M.E., Lazarides, G. and Pallis, C. 2000, Phys. Rev. D 61, 123512. 

Gomez, M.E., Lazarides, G. and PaUis, C. 2000, Phys. Lett. B 487, 313. 

Gomez, M.E. and Vergados, J.D. press, hep-ph/0105115. 

Gondolo P. 2001, Phys. Rev. D 64, 023504. 

Green, A 2001, Phys. Rev. D 63, 103003. 

Gross, M.A.K. Somerville, R.S., Primack, J.R., Holtzman , J. and Klypin, A. A. 1998, Mon. 
Not. R. Astron. Soc. 301, 81. 

Henriksen, R.N. and Widrow L.M., MNRAS, 276, 679. 

Jungman, G. et al.(1996), Phys. Rep. 267, 195. 

Kosmas, T.S. and Vergados, J.D. 1997, Phys. Rev. D 55, 1752; from Yadernaya Fisika, Vol. 
61, No 7, 1166 1998. 

Merritt, D. 1985, A J 90, 1027. 

Perlmutter, S. et al 1999, ApJ 517, 565; astro-ph/9812133. 

Perlmutter, S. et al 1997, ApJ 1997 483, 565; astro-ph/9608192. 

Perlmutter, S., Turner, M.S. and White, M. 1999, Phys. Rev. Let. 83, 670. 

J.R. Primack, astro-ph/0205391. 

Ressell, M.T. 1993, et al., Phys. Rev. D 48, 5519. 

Riess, A.G. et al 1998, AJ 116, 1009. 



-38- 



Sikivie, P. 1999, Phys. Rev. D 60, 063501. 

Sikivie, P. 1998, Phys. Lett. 432, 139; astro-ph/9705038. 

Smoot, G.F. et al. 1992, COBE data, ApJ 396, LI. 

Somerville, R.S., Primack,J R. and Fabers, S.M. press, Mon. Not. R. Astron. Soc. in press; 
astro-ph/9806228. 

Spira, M. et al 1995, Nucl. Phys. B453, 17. 

Turner, M.S. 1990, Phys. Rep. 333-334, 619; astro-ph/9904051. 

Ulho, P., Kamionkowski, M. 2001, JHEP 0103, 049. 

Vergados, J. D. (1990), Supersymmetric Dark Matter Detection- The Directional Rate and 
the Modulation Effect, hep-ph/0010151; 

Vergados, J.D. 1996, J. of Phys. G 22, 253. 

Vergados, J.D. 1998, Phys. Rev. D 58, 103001-1. 

Vergados, J.D. 1999, Phys. Rev. Lett 83, 3597. 

Vergados, J.D. 2000, Phys. Rev. D 62, 023519; astro-ph/0109296. 

Vergados, J.D. 2001, Phys. Rev. D 63, 06351. 

Vergados, J.D. 2001, SUSY Dark Matter in Universe- Theoretical Direct Detection Rates, 
Proc. NANP-01, International Conference on Non Accelerator New Physics, Dubna, 
Russia, June 19-23, 2001, Editors V. Bednyakov and S. Kovalenko; hep-ph/0201014. 

Vergados, J.D. 2002, to be published; astro-ph/0201269. 

Widrow, L.M, astro-ph/0003302. 



This preprint was prepared with the AAS IM^jX macros v5.0. 



-39- 



Table 1: The various parameters describing the dark matter distribution for = 20 (for the 
definitions see text). 
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A 


a 


hy 






0.00 


0.047 


0.000 


4.61 


0.000 




-0.16 


0.064 


-0.010 


3.97 


-0.021 




-0.24 


0.084 


-0.020 


3.50 


-0.038 




-0.30 


0.134 


-0.040 


2.79 


-0.080 


A = 1/2 


1.45 


0.021 


0.031 


6.78 


0.062 




3.34 


0.015 


0.050 


8.02 


0.100 




6.95 


0.011 


0.075 


9.40 


0.150 




12.0 


0.008 


0.100 


10.61 


0.200 




0.00 


0.062 


0.000 


4.05 


0.000 




-0.13 


0.088 


-0.011 


3.41 


-0.022 




-0.20 


0.102 


-0.020 


3.18 


-0.040 




-0.27 


0.171 


-0.046 


2.50 


-0.089 


A = 1 


0.70 


0.037 


0.026 


5.21 


0.051 




2.00 


0.026 


0.052 


6.13 


0.103 




4.20 


0.018 


0.078 


7.30 


0.155 




7.50 


0.013 


0.099 


8.50 


0.201 




0.00 


0.025 


0.000 


6.37 


0.000 




-0.15 


0.032 


-0.005 


5.57 


-0.010 




-0.25 


0.038 


-0.008 


5.15 


-0.015 




-0.40 


0.050 


-0.012 


4.51 


-0.025 


A = -1/2 


2.70 


0.010 


0.025 


10.26 


0.050 




9.20 


0.005 


0.050 


13.22 


0.100 




19.4 


0.004 


0.075 


15.53 


0.150 




34.0 


0.003 


0.101 


17.53 


0.200 
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Fig. 19. — The function T(u) for dark matter in the Eddington theory for the coherent 
mechanism in the case of the target ^^^J. The graphs have been labeled as in Fig. 18. For 
illustration purposes an LSP mass of lOOGeV has been chosen. 



